#Enter data
y <- c(69.3,75.3,75.9,79.3,67.4,74.9,74.4,78.8,73.4,77.5,78.8,81.6,72.3,78.8,78.1,81.5)
#Enter data matrix
Pressure=rep(c(-1,1),8)
Placement=rep(rep(c(-1,1),each=2),4)
Angle=rep(rep(c(-1,1),each=4),2)
Speed=rep(c(-1,1),each=8)
PrxPl=Pressure*Placement
PrxAn=Pressure*Angle
PrxSp=Pressure*Speed
PlxAn=Placement*Angle
PlxSp=Placement*Speed
AnxSp=Angle*Speed
X=cbind(Pressure,Placement,Angle,Speed,PrxPl,PrxAn,PrxSp,PlxAn,PlxSp,AnxSp)
#Fit the model
#Be sure to load BsMD first
library("BsMD")
#Look up BsProb in help to identify the arguments below
ViolinEx=BsProb(X, y, blk=0, mFac=10, mInt = 2, p = 0.25, g = c(2,2), ng = 1, nMod = 10)
summary(ViolinEx)
plot(ViolinEx)